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Abstract 

We discuss the spatial properties of quantum radiation emitted by a multipole 
transition in a single atom. The qualitative difference between the represen- 
tations of plane and spherical waves of photons is examined. In particular, 
the spatial inhomogeneity of the zero-point oscillations of multipole field is 
shown. We show that the vacuum noise of polarization is concentrated in 
a certain vicinity of atoms where it strongly exceeds the level predicted by 
the representation of the plane waves. A new general polarization matrix is 
proposed. It is shown that the polarization and its vacuum noise strongly 
depend on the distance from the source. 
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I. INTRODUCTION 



It is well known that the Jaynes-Cummings model M plays an important role in investi- 
gation of interaction between atoms and quantum radiation field (e.g., see The point 
is that the model describes fairly well the physical processes in the system and, at the same 
time, allows an exact solution. 

In the usual formulation of the Jaynes-Cummings model [jlj-g), the atom is considered 
as though it consists of two or very few non-degenerated levels. In fact, the radiative 
transitions in real atoms occur between the states with given angular momentum j > 1 and 
its projection m = —j, ■ ■ ■ , j (e.g., see This means that even in the case of only two 
levels, the degeneration with respect to the quantum number m taking (2j + 1) different 
values should be taken into account. The simplest example is provided by a dipole transition 
between the states \j = l,m = 0, ±1) and \ f — 0,m' — 0) when the excited atomic state is 
a triply degenerate one (see figure 1). 

Let us stress one more important difference. The radiation field in conventional Jaynes- 
Cummings model is represented by the plane waves of photons with given linear momentum 
and polarization. At the same time, the multipole transitions in real atoms emit the mul- 
tipole photons represented by the quantized spherical waves with given angular momentum 
and parity J7||| . Although there is no principle difference between the plane and spherical 
waves within the classical domain since both represent the complete orthogonal sets of so- 
lutions of the homogeneous wave equation |J and can be re-expanded with respect to each 
other, the quantum counterparts of these two representations are non-equivalent because 
they describe the physical quantities (the linear and angular momenta respectively) which 
cannot be measured at once. 

Moreover, the multipole field is characterized by more quantum degrees of freedom than 
the plane waves of photons. In fact, the monochromatic pure j-pole multipole radiation of 
a given type (either electric or magnetic) is specified by the (2j + 1) different values of the 
quantum number m. Since j > 1, the total number of degrees of freedom here is not less 
than three. At the same time, the monochromatic plane waves of photons are described 
by only two different polarizations. In particular, the increase of the number of degrees of 



freedom can lead to an increase of the zero-point oscillations [1C 



The multipole generalization of the Jaynes-Cummings model has been discussed in 



11,12]. Let us stress that similar models have been considered in different problems of 



interaction of quantum light with matter (e.g., see jl3Hl6|l and references therein). 

The main objective of this paper is to examine the quantum properties of light emitted 
by a dipole atom at any distance from the source, depending on the boundary conditions. 
The paper is arranged as follows. In section 2 we review the properties of quantum multipole 
field in comparison with these of plane waves of photons. In section 3 we briefly discuss the 
multipole generalization of the Jaynes-Cummings model. Then, in section 4 we consider 
the polarization of multipole radiation and introduce a novel general polarization matrix. 
This object permits us to take into account the spatial anisotrophy of both the electric and 
magnetic fields at once. In section 5 we examine the spatial properties of multipole photons 
emitted by an atom in an ideal cavity as well as in empty space. In particular, we show that 
the polarization properties of quantum multipole radiation changes with distance from the 
atom. In section 6 we briefly discuss the obtained results. 
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II. QUANTUM MULTIPOLE FIELD 



Following [f?|,jS|jr7l , we list below some important formulas describing the quantum mul- 
tipole field. It is usually considered in the so-called helicity basis HT7 | 



X± = =F ^ , Xo = e z . (1) 

It is clear that {x^} formally coincides with the three eigenstates of spin 1 of a photon. Since 
the polarization is defined to be the spin state of photons ||18|| , one can choose to interpret 
X± as the unit vector of circular polarization with either positive or negative helicity, while 
Xo gives the linear polarization in the z-direction. To within the sign at x± the helicity basis 
(1) coincides with the so-called polarization basis usually used in optics [19|]. In the basis 
(1), the positive-frequency part of the operator vector potential of multipole field can be 
expanded as follows |7|Jl7|| 

1 oo j 

Mrl = E E E E (-lYx-»V Xkjm ,(7)a xkjm , (2) 

k n=-lj=lm=-j 

where a... is the photon annihilation operator which obeys the following commutation rela- 
tions 

Here k is the wave vector, A = E,M denotes the type of radiation (parity), index j > 1 
gives the angular momentum, index m = —j, ■ ■ ■ , j, and a\kj m is the annihilation operator of 
the corresponding photon. The mode functions V...(r) in (2) are represented in the following 

way 



V E kjmn = lEkj[\fjfj+i(kr)(l,j + 1, fj,,m - fJ,\jm)Y :j+lim _ ll 

- yj + lfj-i(kr)(l,j - l,/i,m-/i|jm)F j _ lim _ M (^,0)], 
V M kj mf i = lMkjfj(kr)(l,j, fi,m- fi\jm)Y jtm _^(9, 4>) (3) 

for the electric and magnetic radiation respectively. Here 



27rhc 2irhc 



lEkj — \ , , ,/ . -. v iMkj 



fcV(2j + l)' ,m " J V kV ' 

are the normalization constants, V is the volume of quantization, (■ ■ ■ \jm) denotes the 
Clebsch- Gordon coefficient of vector addition of the spin and orbital parts of the angular 
momentum of a multipole photon, and Yt m is the spherical harmonics. The radial contribu- 
tion into the mode function (3) depends on the boundary conditions. In the standard case 
of quantization in terms of standing spherical waves in a spherical cavity |7,I7|, we have 



fi(kr) = j e (kr) = JJLj t+l/2 (kr), (4) 
where je(x) is the spherical Bessel function. 



3 



The positive frequency parts of the operator field strengths obey the following relations 

EEkjm =i E f t (~^) fl X-^Ekjmfi(r) a Ekjm, 
fj.,k,j,m 

EMkjm — i E ^)^X-/jVMkjmti('^} a Mkjm, 
ti,k,j,m 

BEkjm = ~ i E k( — l) fl X~^Mkjm^( : f f )0'Ekjm, 
fi,kj,m 

B M kjm = i E K-^YX~^yEkjrnpi{f)a Mkjm . (5) 

It can be easily seen from (2) and (5) that the electric multipole field always has longitudinal 
component of the electric field strength in addition to the two transversal components, 
while it is completely transversal with respect to magnetic induction. At the same time, 
the magnetic multipole field has all three components of magnetic induction and only two 
transversal components of the electric field strength. 

The position dependence of the mode functions (3) is not an unusual fact. In reality, the 
mode functions of the plane waves also depend on position: 



A {plane) = Ey ^ E(-l) M X-^ ifc> a fc ,- (6) 

Here we choose the basis (1) with xo — k/k and a... denotes the photon annihilation operator, 
corresponding to the states with given linear momentum (direction of propagation) and 
transversal polarization with either helicity. The third projection of the photon spin is 
forbidden in this case [(18] . 



To complete the discussion of the quantum multipole field, consider now the zero-point 
or vacuum oscillations. Since the free-field Hamiltonian has the form [|7|J§[17| 

H = E huj k E (alkjm a Xkjm + 1/2), 

k 

for the energy of the vacuum state in the whole volume of quantization we get 

H vac = E^ fc E 1/2 = E^ (E(2i + i) ) ■ (7) 

k Aj',m k \ j J 

For comparison, we show here well known expressions valid in the case of plane waves of 
photons (6): 

^=E^(a^ + l/2), 

k,fi 

H ( P lan e) =Y^kY. I/ 2 = E (8) 
k £t=±l k 

Due to the definition of k, both expressions (7) and (8) give an infinite energy and, at 
first sight, cannot be compared with each other. In fact, this infinity is inessential because 
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of the following reason. The contribution of zero-point oscillations can be observed only 
via measurement which implies an averaging of physical quantities over a finite "volume of 
detection" and exposition time of detector • Such an averaging plays a part of filtration 

leading to a selection of a certain finite transmission frequency band. It is therefore seen 
from (7) and (8) that, even if the filtration process leads to selection of the electric dipole 
photons only, the ratio of contributions of the zero-point oscillations of the multipole field 
and plane waves is equal to 3/2. From the physical point of view, this result is caused by 
the more number of quantum degrees of freedom in the case of multipole photons |10 . 



Much more interesting and important result can be obtained from consideration of the 
spatial properties of the zero-point oscillations. The energy density of the field is represented 
as follows 

107T 107T L > 

Then, the zero-point contribution is 

(0\m\0) = ^ {(0\E ■ E + \0) + (0\B- B+\0}} . (9) 

Consider first the case of plane waves. In view of (6) and symmetry relations E x = 
B y , E y = -B x , we get 

m iplane \m = ^T. h ^. (io) 

k 

Thus, the zero-point oscillations of the energy density of the plane waves of photons are 
homogeneous in the space of quantization in spite of the position dependence of the mode 
functions in (6). 

In turn, employing the relations (5), we get 

(o|#(r)|o) = ^-£* 2 E l*W*OI 2 - (ii) 



8tt , 

k A,j,m,[i 



In view of (3), we obtain 



(0\H(r)\0) = -|t^H E [^-^lv^/i+i( A;r )(l'^ + l'/ i ' m -/ i b' m ) y i+i« 

-\Jj + 1/j_i(At)(1, j - l,/i,m -/i|jm)r i _ 1>m _^(6',0)| 2 

+\f j (kr)(lJ,^m-v\jm}Y j:m ^(e,(f))\ 2 ]. (12) 



Taking into account that 



it is straightforward to show that the zero-point energy density of the multipole field is 
independent of the angular variables 9 and <p, while depends on the distance from the origin 
(the singular point corresponding to the source location): 
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(0\H(r)\0) = (0\H(r)\0). 

Thus, the density of the zero-point oscillations of multipole field has the spherical symmetry 
with respect to the singular point (atom), while manifests the radial dependence. Taking 



into account the radial dependence (4) and making use of the following formula [22 

it is easy to prove that the multipole zero-point oscillations (12) vanish at far distance. Thus, 
the vacuum fluctuations of the multipole field are concentrated near the atom. 

Further employing the properties of spherical Bessel functions shows that the principal 
contribution into (12) at kr — > comes from the term with jo(kr), corresponding to the 
electric dipole radiation. The radial dependence of (12) at fixed k and j — 1 is shown in 
figure 2. It is seen that the zero-point oscillations of the multipole field are concentrated 
in some vicinity of the singular point where they strongly exceed the level (10) predicted 
within the framework of the representation of plane waves of photons. From the figure 2, 
the radius of the region of concentration of the zero-point oscillations can be estimated as 
follows: 

2 Afc Afc . . 

ro ~fc = T~y (13) 

where is the wavelength. Thus, the atom condenses the zero-point oscillations in the 
near and intermediate zones. Let us stress that in the number of modern experiments on 
engineered entanglement in the systems of trapped Rydberg atoms, the interatomic distances 
are of the same order ||23|| . 

The above result (13) should be considered as an estimation from below because it 
corresponds to the first term in (12). Successive taking into account of further terms leads 
to a certain shift of r into the intermediate zone. 



III. MULTIPOLE JAYNES-CUMMINGS MODEL 

The main aim of this section is to emphasize that the interaction between the photons 
and electrons on atomic sub-levels with the same j and different m is specified by a certain 
coupling constant independent of m. 

Consider a two-level atom with the electric dipole transition j = 1 — > j' = 0. The cou- 
pling constant of the atom-field interaction can be found by calculating the matrix element 

g = — (0,0\p- A + A-p\l,m) = ik o (0, Old- All, m), (14) 

2m e c 

obtained from the expression 
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describing the interaction between the atomic electron with linear momentum p, charge e, 
and mass m e and radiation field specified by the vector potential A. Here d = er is the 
dipole moment of the atomic transition with the resonance frequency u = ck . Assuming 
the central symmetry of atomic field and taking into account the fact that the spin state of 
an atom does not change under the electric dipole transition, we can represent the atomic 
states in (14) as follows 

((r)\l,m)=n e (kr)Y lm (0,<p), ((f)\0, 0) = K g (kr)Y 00 (9 , 0), 

where lZ e ,g is the radial part of the atomic wave function of either excited (e) or ground (g) 
state. 

Expanding the dipole moment d over the helicity basis (1), substituting (2), and carrying 
out the calculations of integrals in (14) over a small volume occupied by the atom, for the 
coupling constant (14) we get 

9=^D, (15) 

V KC 



where D is the effective dipole factor which, by construction, is independent of the quan- 
tum number m. Taking into account the properties of the Clebsch- Gordon coefficients and 
spherical harmonics, for the position-dependent mode function in (2) we get 

lirn Vskimn ~ <W- 

This means that the electric dipole transition |1, m) — > |0, 0) at any given m creates a photon 
with helicity (polarization) fi = m — 0, ±1. 

Finally, the Jaynes-Cummings Hamiltonian of the electric dipole transition in the 
rotating-wave approximation [^4j takes the form [|TT|,|r2] 

fT H = H Q + H int , 

i 

Ho = ^ {uja^cim + uj R mm }, 

m=—l 
1 

Hint = g X! {Rrngdm + a m Rgm} ■ (16) 
m=— 1 

To simplify the notation, hereafter we omit insignificant indices. Here the atomic operators 
are defined as usually ||23] in terms of the projections on the atomic states: 

R mg = |l,m)(0,0|, R mm i = |l,m)(l,m'|. 

The Hamiltonian (16) describes the creation and absorption of the single cavity-mode pho- 
tons at the atom location. Everywhere in the surrounding space, we have to take into account 
the spatial dependence of the radiation field described by the vector potential (2). Similar 
model can be constructed in the case of magnetic dipole radiation as well as in the case of 
other high-order atomic multipoles. 
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IV. POLARIZATION OF MULTIPOLE RADIATION 



It is known that the polarization defines the direction of oscillations of the field strengths. 
Within the classical picture based on the consideration of plane waves, the polarization is 
defined to be the measure of transversal anisotrophy of the electric field strength [|L9[| . In turn, 
the quantum mechanics interprets the polarization as given spin state of photons ||18|| . In the 
usual approach, the quantitative description of polarization is based either on the Hermitian 
Polarization matrix or on the equivalent set of real Stokes parameters. In the standard case 
of plane waves, we get the (2 x 2) polarization matrix and four Stokes parameters ITS, while 
the description of multipole radiation requires for the (3 x 3) polarization matrix and nine 
Stokes parameters |2^J2^]. Moreover, the electric- and magnetic-type radiation fields are 
usually described in terms of different polarization matrices, taking into account the spatial 
anisotrophy [E7|. 

Here we construct a more general novel object, describing in a unique way the polarization 
properties of multipole radiation of either type both classical and quantum as well as these of 
the plane waves and other forms of electromagnetic radiation (e.g., of the cylindrical waves). 
In general, the field can be described in terms of the field-strength tensor which can be 
2§ 



chosen as follows 



/ 





E x 


Ey 


E z \ 




-E x 





-B z 


By 




~Ey 


B z 





-B x 


V 


-E z 


-By 


B x 


o / 



(17) 



It seems to be tempting to introduce the general quantitative description of polarization us- 
ing (17). Since the polarization is specified by the intensities of different spatial components 
of the radiation field and by the phase differences between these components |19|], it should 
be described in terms of bilinear forms in the field strengths. The simplest bilinear form in 
the field-strength tensor is 



R = F+F, 



(18) 



which differs from the energy-momentum tensor by a scalar. In some sense it is similar to 
the Ricci tensor considered in the general relativity |2"9" |. It is easily seen that (18) has the 
following structure 



R 




where We = E + ■ E is a scalar, S, apart from an unimportant factor, coincides with the 
Poynting vector, and P is the Hermitian (3 x 3) matrix of the form 



P = P f 



Pf. 



Here 



/ E+E x E+E y E+E z 
E+E x E+Ey E+E z 
V EtE x EtEy EtE z 



(19) 



(20) 



S 



and 



/ B+By + B+ Z B Z 



B 



V 



-BtBy 

-BtB, 



-B+B x -B+B x 
B+B x + B+B z -BtB, 
-B+B z B+ X B X 



B } B y 



(21) 



We note here that (20) has been proposed in [EH] in order to describe the spatial anisotrophy 
of the electric dipole radiation, while (21) is similar to the object has been discussed in 

2is 

We choose to interpret (19) as the general polarization matrix, while the terms (20) and 
(21) give the electric and magnetic field contributions respectively. 

To justify this statement, consider first the case of plane waves propagating in the z- 



direction when E z 
form 



and B x = —E y , B y = E x . Then, the matrix(20) takes the 



j(piarae) 



/ B+E x E+E y 
E+E x E+Ey 
V 



(22) 



It is seen that the non-zero submatrix in (22) coincides with conventional (2 x 2) polarization 
matrix of plane waves [19,21]. In turn, (21) takes the form 



Pf. 



( E+E x E+ x E y 
E+E x E+Ey 



(23) 



V o 



E + ■ E 



where corresponding (2 x 2) submatrix in the top left corner coincides with (22). Thus, the 
general polarization matrix (19) describes the polarization of plane waves adequately. 

Consider now the multipole radiation. In the case of electric-type radiation when B z = 
everywhere, the matrix P% has the general form (20), while the magnetic polarization matrix 
(21) is reduced to 



P, 



B 



( B+B y -B+B x 
-B+B y B+B x 



V 













B+ ■ B 



(24) 



The polarization of magnetic-type radiation is described by (20) with E z = which coincides 
with (22) and by the general form (21) with B z ^ 0. In this case, (21) coincides, within the 
transposition of lines and columns, with the polarization matrix considered in j27 . 



It is natural that the general polarization matrix (19) reflects the three-dimensional 
structure of the radiation field. The diagonal terms in (20) and (21) give the radiation 
intensities. Their angular and radial dependence corresponds to the radiation patterns of 
the multipole field. The off-diagonal terms give the phase information as in the case of plane 
waves [0. In contrast to the standard case of plane waves, there are the two independent 
phase differences Ajj = arg E± — arg Ej instead of only one phase difference because 



A 



A, 



0. 
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Since E(r) ■ B(r) = at any point, the magnetic part (21) of the general polarization matrix 
(19) contains the same phase differences as (20). 

Similar expressions for the polarization matrix (19) can also be obtained in the helicity 
basis (1). For example, the matrix (20) takes the form 



( E+E+ -E+E E+E_ \ 

-E+E + E+E -E+E_ 
\ E+E+ -E±E E+E_ ) 



(25) 



The quantum counterpart of (19) can be obtained by formal substitution of the operators 
instead of the classical field strengths (compare with ||30|| ). Averaging of the correspond- 
ing operator matrix over a given state of the radiation field gives the polarization matrix. 
By construction, the operator matrices (19)-(24) correspond to the normal ordering in the 
creation and annihilation operators: 

P = P( a + a). 

In addition, one can define the anti-normal polarization matrix 

p(an) = p( fla +) 

by a simple change of order of product of the field strengths in all elements of the matrices 
(19)-(24). It is then clear that the matrix 

p(an) _ p = P ([ a ,a + ]) = (0|P^|0) EE P vac (26) 

determines the zero-point (vacuum) contribution into the polarization. Following the ideas 
and results of section 2, it is a straightforward matter to show that the vacuum polarization 
of plane waves of photons is uniform in the space, while the multipole vacuum polarization 
concentrates near the atom and exceeds the level predicted by the representation of plane 
waves. 



V. POLARIZATION OF A SINGLE- ATOM RADIATION 

It was shown in section 3 that an atomic electric dipole transition with given m emits 
the photon with given polarization fi = m. We now show that the polarization changes with 
the distance from the source. In other words, the polarization is not a global property of 
the field, while changes from point to point, at least in the case of multipole radiation. 

Assume, for example, that the atom emits the electric dipole photon with m = +1, i.e. 
circularly polarized with positive helicity. Consider the polar direction (6 = 0), correspond- 
ing to the maximum of the radiation pattern in this case |J. Then, the matrix (25) averaged 
over the photon state |1 + ) takes the form 

[\h{kr) - j (kr)} 2 \ 
Pr = — I 

0/ 
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Thus, there is only one polarization // = +1 in the polar direction. In the less probable case 
of the equatorial direction (9 = ir/2), from (25) we get 



huj I llh(kr) + jo(kr)}* =fj 2 (kr)[\j 2 (kr) + j (kr)]e** \ 



3V 





\fj2(kr)[lj 2 (kr)+ Jo (kr)]e- 2 ^ U^( kr )f ) 



(27) 



so that there are the two circularly polarized components with opposite helicities. Compar- 
ison of intensities of the two components shows that the positive helicity dominates at short 
distances (kr < 3), while both components contribute equally at far distances (kr ^> 1) 
(see figure 3). It is clear that any deviation from the polar direction leads to creation of 
polarizations additional to fi = +1. Thus, the polarization o radiation under consideration 
strongly depends on the direction and distance from the source. Similar picture can be 
obtained for polarization of photons with m — — 1 and m — 0. 

The above results were obtained in the case of standing waves in an ideal spherical cavity 
when the radial dependence of the mode functions (3) is specified by equation (4). In this 
case, the radiation field is subjected to the Rabi oscillations which can be described through 
the use of the steady-state time dependent wave function for the system with Hamiltonian 
(16): 

\m) = e- tHt \e m , 0) = \ £ e-^*(|e m ; 0) + £\g; l m », (28) 

A e=±i 

where we choose the initial state as the vacuum state of the cavity field and excited state of 
the atom is specified by given m. Then, the elements of the polarization matrix (19) obtained 
by averaging or corresponding operator matrix over the state (28) should be multiplied by an 
additional factor of (1 — cos 2gt), describing the steady-state time dependence of polarization. 

Consider now the radiation by a dipole atom in empty space. Then, the Hamiltonian 
(16) should be generalized as follows 

i 

Hp — ^ {y^ j UkQ>tm. a km + UpRmm}, 
m=— 1 k 
1 

Hint = X! 9k 51 {Rmg a k in + at m R gm }, (29) 
k m=— 1 

to take into account the /c-dependence of the radiation field. Let us again choose the initial 
state as the vacuum state of photons and excited atomic state with given m 

hM = |em)®[®|0 fc )]- (30) 
k 

It is then clear that the radiation field is represented by the outgoing spherical waves of 
photons which stipulates the choice of the radial dependence in (3) in terms of the spherical 
Hankel function of the first kind 

fe(kr) = h { l\kr) = Je (kr) + i{-lf +1 j ^{kr) (31) 
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instead of (4) ||31|| . This choice assumes that the atom occupies a small but finite spherical 
volume of radius r a at the origin to avoid the divergence at kr — > 0. 

The elements of the polarization matrix (25) of the electric dipole radiation in the equa- 
torial direction take the form 

EXE + = ^{[T + {kr)f + [T.{kr)f}, 

ElE^ = ^{\j 2 (kr)r + \j_ 3 (kr)r}, 

EXE- = ^{[j 2 (kr)T + (kr)+j. 3 (kr)T4kr)] 2 

+ [j 2 (kr)T4kr) - U(kr)r + (kr)] 2 } l / 2 e^, (32) 



where 



F+(kr) = ^j 2 (kr) + j (kr) 
T-(kr) = ^j_ 3 (At) + j_i(fcr) 



and 



ip = 20 + tan 1 



V-{kr)j 2 (kr) - T + (kr)j- 3 (kr) 
T + (kr)j 2 (kr) + Y -{kr)j- 3 {kr) 



All elements containing E are equal to zero. Unlike the case of radiation in an ideal cavity 
(27), the phase difference <p between the components with opposite helicity depends here on 
the distance from the source. Other cases can be examined in the same way. 

For the polarization matrix (25) in the polar direction, we again get only one non-zero 
element 

E + + E + = ^ {[E + (kr]} 2 + £-(kr)] 2 } , 

where 

Z+(kr) = ^j 2 {kr) - j (kr) 
E_(Ar) = -j- 3 (kr) -j_i(Ar) 



The time evolution can be described with the aid of approach proposed in |32j| . Then, the 
matrix elements of the polarization matrix should be multiplied by the following factor 

1 - 2e-" i/2 cos[(cu fc - uj - Au)t] + e" 7 *, 

where 

7T V 

* k ^0 ^k 

and V denotes the principle value of corresponding integral. Again, any deviation from the 
polar direction changes the picture of polarization at intermediate and far distances with 
respect to that at the atom location. 
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VI. CONCLUSION 



Let us briefly discuss the obtained results. In this paper we have concentrated on the 
description of the spatial properties of multipole radiation by a single atom. The consid- 
eration is based on sequential use of the representation of multipole photons corresponding 
to the radiation of real atoms. It is shown that the zero-point oscillations of the multipole 
field are concentrated near the atom, while exceed the level predicted by the model of plane 
waves of photons everywhere. In a certain neighborhood of the atom, the effect is strong 
enough. Although the effect can be observed at the relatively short distances, it seems to be 
important for the near- and intermediate-field quantum optics as well as for the experiments 
with trapped Rydberg atoms when the typical interatomic distances are of the same order 
P5fl . In particular, it can be important for engineered entanglement in the system of two 
atoms in a cavity proposed in |H| , for experiments with single-atom laser as well as for the 
estimation of Casimir effect in atomic systems. 

Let us note in this connection that possible influence of an atom on the electromagnetic 
vacuum state in the absence of radiation has been discussed in quantum electrodynamics 
for a long time (e.g., see B3J3H]). The new element here is the spatial inhomogeneity of the 



vacuum noise. Unlike the effects discussed in [j34[ , j3~5| , the specific distance dependence of the 



zero-point oscillations in the presence of atom has the geometrical nature and is independent 
of the atom-field interaction. 

To describe the polarization of multipole field, we proposed in section 4 a new definition 
of the polarization matrix based on the bilinear form in the field-strength tensor. The 
generalized polarization matrix (19) is additive with respect to the contributions coming 
from the electric field and magnetic induction. It reflects the three-dimensional nature of 
polarization connected with the three possible states of spin of photon. In special case of 
plane waves, when the third spin state is forbidden, it reduces to the conventional (2 x 2) 
Hermitian polarization matrix. In the case of multipole radiation, it combines together the 
objects considered earlier in poVpTj . 

By construction, the generalized operator polarization matrix is a local object in spite of 
the global nature of the photon operators of creation and annihilation. The spatial properties 
of polarization are caused by the mode functions and changes with distance and direction 
from the source. For example, the electric dipole radiation from the excited atomic level 
with m — +1, corresponding to the creation of a circularly polarized photon with positive 
helicity, is transformed, at far distances, into the radiation with both helicities. In the case 
of radiation in empty space, the phase difference between the modes with opposite helicities 
is also a function of distance and direction from the source. 

Both the distance and direction dependence of the polarization seems to be very im- 
portant. In fact, any real measurement of intensity assumes the finite aperture of a de- 
tecting device |21]| , |2"T| ] . This means that, in the case of radiation by the atomic transition 
j = l,m = +1) — > \j' = 0, ml — 0) considered in section 5, the precision of measurement 
of the polarization /x = +1 in the polar direction would be influenced by the zero-point 
oscillations of all three polarizations. 

Let us stress that locality of polarization discussed in this paper can also be interpreted 
in terms of the photon localization. It is well known that, while the photon operators 
of creation and annihilation are defined in the whole space at once, the notion of photon 
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localization is not deprived of physical meaning [pl|,j34 
a quite certain example of the photon localization 
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The photodetection process provides 
. There are also known attempts to 
interpret the photon localization as the specific fall-off of the photon energy density P7I . 
It is also clear that emission and absorption of radiation by atoms can also be interpreted 
as the photon localization ||38|| . Our results based on consideration of the representation of 
multipole photons shows the rapid fall-off of the density of zero-point oscillations and of the 
vacuum noise of polarization with the distance from atom as well as the spatial dependence 
of a certain characteristics of the radiation field (e.g., polarization). 

One of the authors (A.S.Sh.) would like to thank Professor J.H. Eberly, Prof. A. A. 
Klyachko, Professor V.I. Rupasov, and Professor A. Vourdas for many fruitful discussions. 
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FIGURES 

FIG. 1. Energy diagram of triple degenerated excited and ground states of a dipole transition 
j = l~j' = 0. 

FIG. 2. Zero-point oscillations of the energy density versus dimensionless distance kr at fixed 
k and j = 1. 

FIG. 3. Distance dependence of intensity of the multipole radiation generated by the atomic 
transition \j = l,m = +1) — > \f = 0, m! = 0) in equatorial direction. 
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